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TECHNICAL NOTE 3058 


TRANSIENT TEMPERATURE DISTRIBUTIONS IN SIMPLE 
CONDUCTING BODIES STEADILY HEATED THROUGH 
A LAMINAR BOUNDARY LAYER 
By Hermon M. Parker 


SUMMARY 


An analysis is made of the transient heat-conduction effects in 
three simple semi -infinite bodies: the flat insulated plate, the conical 

shell, and the slender solid cone. The bodies are assumed to have con- 
stant initial temperatures and, at zero time, to begin to move at a con- 
stant speed and zero angle of attack through a homogeneous atmosphere. 

The heat input is taken as that through a laminar boundary layer. Radi- 
ation heat transfers and transverse temperature gradients are assumed to 
be zero. 

The appropriate heat -conduction equations are solved by an iteration 
method, the zeroeth-order terms describing the situation in the limit of 
small time. The method is presented and the solutions are calculated to 
three orders which are sufficient to give reasonably accurate results 
when the forward edge has attained one-half the total temperature rise 
(nose half -rise time). Flight Mach number and air properties occur as 
parameters in the result. Approximate expressions for the extent of the 
conduction region and nose half -rise times as functions of the param- 
eters of the problem are presented. 


INTRODUCTION 


One of the major problems arising in supersonic and hypersonic 
flight of aircraft is that of aerodynamic heating. It is apparent that 
the maximum feasible speed for a given aircraft at a given altitude is 
l im ited by this effect. This problem has been studied from many varying 
points of view. Usually (refs. 1 to 3 ) these investigations were con- 
cerned with determining the temperature attained at the surface of a 
given body in steady flight or during a given flight history of variable 
speed and perhaps variable altitude. Little attention has been devoted 
to the effect of heat conduction within the body on the transient temper- 
ature distribution of its surface. Kaye (ref. 2) considered the aero- 
dynamic heating of an infinite-aspect-ratio, 8-percent-thick, solid-steel. 
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symmetrical double -wedge wing of 5 • 5-foot chord in flight in a homogeneous 
atmosphere and accelerating at a constant rate from Mach number 1.4 to 
Mach n umb er 6.0. Using a "coarse network" of integration points, he 
numerically integrated the appropriate two-dimensional heat-conduction 
equations, using the best available experimental values of the heat- 
transfer coefficient (ref. 4) and neglecting radiation transfers. He 
concluded that, except near the leading and trailing edges, the chord- 
wise conduction is relatively unimportant but that the transverse temper- 
ature gradients, especially at midchord, are of considerable importance 
in estimating the transient surface temperature distribution. 

For actual aircraft in high-speed flight, the maximum rate of rise 
of temperature occurs at the leading edge or nose. Under the assumptions 
of zero internal, conduction, of no heat transfer except that through the 
boundary layer, and of steady flight beginning at zero time in a homo- 
geneous atmosphere, the surface temperature approaches the local adia- 
batic wall temperature at a rate which is theoretically infinite at zero 
Heynolds number and decreases with increasing Reynolds number. With 
internal conduction the surface approaches an equilibrium temperature 
distribution at a rate which is finite everywhere and is appreciably 
different from the zero- conduction rate at the smaller Reynolds numbers. 
Nonweiler (ref. 3) has calculated the equilibrium temperatures near the 
leading edge of the flat plate in steady flight by taking into account 
conduction in the plate, heat input through a laminar boundary layer, 
and loss of heat by radiation with the assumption that the difference 
between the local equilibrium temperature and the adiabatic wall temper- 
ature is large and constant. 

The temperature lag due to heat capacity and heat conduction permits 
aircraft to exceed for short times the speed for which the corresponding 
equilibrium temperature is the maximum permissible temperature of the 
body. Further, transient thermal stresses axe important problems in the 
structural design of high-speed aircraft. Thus it is important to deter- 
mine quantitatively the effect of heat conduction on the transient temper- 
atures near the leading edge or nose of as general a class of bodies and 
flight conditions as possible. 

The problem considered herein is that of determining the transient 
effect of heat conduction parallel to the airstream and within certain 
simple infinitely long bodies in steady flight at zero angle of attack 
in a uniform atmosphere. The bodies considered are the insulated flat 
plate, the conical shell, and the slender solid cone. Heat transfers 
due to radiation are neglected. Because the effect of heat conduction 
is Important only near the leading edge or nose, the usual expressions 
for heat transfer through an incompressible laminar boundary layer are 
used. Corrections to the incompressible heat transfer which do not 
alter the temperature and Reynolds number dependence may be included in 
the analysis. Because the lami nar boundary-layer equations do not apply 
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exactly at the leading edge or nose, the results of the analysis are 
applicable for Reynolds number (based on distances from the leading edge) 
greater than the order of 10. Temperatures within the body perpendicular 
to the stream direction are assumed uniform. 
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conduction parameter, [- — — ) — , 


— , sec - V^ 


,1/3 


0.33kgU pr- 1 . 

heat-transfer constant for flat plate, , sec _J - 


Vb' 


heat-transfer constant for slender solid cone, 

5 1/3 
0.33k ITPr 


v h 


& 


g, sec 


-1 


PbV 

speed of sound, ft/sec 
specific heat of material of body 


sin 9 

geometry factor for slender solid cone, 

1 - cos 9 

thickness of flat plate or conical shell, ft 

thermal conductivity of air, Btu/(sec)(ft)(°R) 

thermal conductivity of material of body, Btu/(sec)(ft)(°R) 

Mach number just outside boundary layer, U/c 
summation variable 
Prandtl number 

rate of heat transfer through laminar boundary layer on flat 
plate, Btu/(sec)(ft^) 
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rate of heat transfer through 
cone, q.^3, Btu/(sec)(ft2) 


boundary layer on 


„ , Ux Us 

Reynolds number, — or — 

Reynolds number at rear of conduction region at nose half -rise 
time 

distance variable on conical surface, measured from vertex, ft 
temperature of body, °R 
adiabatic wall temperature, °R 

initial temperature of body, °R 
time, sec 

nose half -rise time, sec 

stream speed just outside boundary layer, fps 

stream speed in boundary layer (Blasius flat-plate solution), fps 


distance variable on flat plate, measured from leading edge, ft 
distance perpendicular to wall (Blasius flat-plate solution), ft 
dimensionless time parameter 
dimensionless distance parameter 


dimensionless transient-heat-conduction parameter. 


at " 1 -/ 2 


value of Tj corresponding to rear boundary of heat-conduction 
region 

density of material of body 
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0 half vertex angle of slender solid cone 

v kinematic viscosity 

Primes with a symbol denote differentiation with respect to t). 
Subscripts: 

x at distance x 

s at distance s 

h for body material 

g for gas (air) 


METHOD OF ANALYSIS 
Rate of Heat Transfer 


For the simple semi-infinite bodies to be considered in this analysis, 
the effect of internal heat conduction parallel to the stream direction 
is significant only near the leading edge or nose. Therefore, it is 
natural to assume for the heat transfer those expressions appropriate to 
a laminar boundary layer. The present analysis is applicable to cases 
in which the expression for the heat transfer to the body is of the form 


q. = 


(Constant) 


T, 


aw 


ft 


or 


q. = 


(Constant ) 



where x is the distance from the leading edge of the flat plate and s 
is the distance along the conical surface from the vertex of the conical 
shell or the slender solid cone. Thus, the real assuuption regarding the 
heat input is the temperature and distance dependence. Compressibility 
and nonuniform-surface -temperature effects which do not alter the temper- 
ature and distance dependence may be included in the present analysis by 
suitable changes in the constant of proportionality. For simplicity the 
usual expression for the heat transfer in the incompressible isothermal 
case is assumed. 
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The rate of heat transfer through a unit area of an incompressible 
laminar boundary layer on the flat plate is 

1 - 0.35kg 2 PrV3 (1) 

®x V 

where R x is based on the distance x from, the leading edge. Equa- 
tion (l) may be derived from the first equation on page 626 of refer- 
ence 5 . Hantzsche and Wendt (ref. 6 ) have shown that the laminar 
boundary-layer heat transfer on a cone is \/3 times that on the flat 
plate with the same free-stream conditions just outside the boundary 
layer. Thus it follows that the rate of heat transfer through unit area 
of a laminar boundary layer on the cone is 


= 0.33 \J3 K % Pr 1 /3 J_? 


( 2 ) 


where R s is based on the distance s from the vertex along the surface 

of the cone. Since no satisfactory representative length exists for the 
semi -infinite bodies to be considered in the analysis, the Reynolds num- 
ber is chosen as a suitable dimensionless distance parameter. 


The expressions (l) and (2) are based on the assumption that the 
Blasius flat-plate solution for the laminar boundary layer is correct. 
Inasmuch as it is known that the Blasius solution is not correct at very 
low Reynolds numbers, an estimate of the Reynolds number at which the 
Blasius solution fails was made by the following calculation: 


8^u 

The Blasius solution was used to calculate the term v — -, which 

Sbc 

was neglected in the Navier-Stokes equations, on the basis that it is 

^u 

s m all in comparison with v . With the curve y(R x ) defined as that 

by 2 

^2 v2 

curve under which v — ^ is equal to or larger than v — it was 


8 x c 




found that y(R x ) equals the displacement thickness at R x 5 and, 
at large R^ is approximately proportional inversely to the 
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square root of R x . Thus a reasonable conclusion appears to be that the 
Blasius solution and therefore equations (l) and (2) should not be in 
serious error for Reynolds numbers larger than the order of 10. 

Equations (l) and (2) indicate that, at the lead in g edge or nose, 
the rate of heat transfer becomes infinite, which, of course, is physi- 
cally impossible. Actually, at small Reynolds numbers the heat transfer 
must be limited by the finiteness of the random molecular motion, as 
happens in free -molecule flows. Comparison of the heat transfer in the 
free-molecule regime with equation (l) leads to the conclusion that the 
maximum q occurs at a Reynolds number of the order of 1/10 of the 
square of the Mach number. On the other hand, the failure of the Blasius 
solution at small Reynolds numbers is surely related to other approxi- 

^2 

mations, such as the neglect of the term v — -. In any event, the 

cbc 

effect of internal conduction is important over a range of Reynolds num- 
bers very large compared with 10, so that the correctness of equations (l) 
and (2) for Reynolds numbers less than 10 is not of great importance. 

The integrability of the heat input assures no mathematical difficulties 
and no serious errors even at reasonably small Reynolds numbers. 


The Flat Plate 

The general three-dimensional heat -conduction equation is (ref. 7): 


v* I = £(*» I) + £(*» I) + !(*»> I) + Q(x ’ y ’ z) 


where k^ is the thermal conductivity and Q(x,y,z) is a heat-generation 

function specifying the rate at which heat is created per unit volume at 
the point (x,y, z). The problem of the flat plate heated through a bound- 
ary layer reduces to the one -dimensional time -dependent heat-conduction 
problem when temperature gradients normal to the stream direction are 
assumed to be zero. The heat-generation function becomes the rate of 
heat transfer per unit area q divided by the plate thickness h which 
is assumed to be constant. Thus, the temperature of the flat plate 
heated from one side is given by 


St a 2 T , „ ,,, U Pr 1 ' 5 

‘Vt. - 'S 13 + °- 55k e V T 


cbc £ 


T, 


aw 


R-v 


- T 

V2~ 
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or 


where 


c& _ 2 cFk , 
^ ' 8r x 2 


— — (1 - K ) 

R x 1/2 


a 2 = ^ ^ 

PbS v 2 


h 


0.33 


k g U Pr 1 / 5 

%% V b 



(3) 


Boundary and initial conditions are chosen as follows. At zero 
time the plate has a uniform temperature T ± so that 


Kt=o =0 PO 

and at large values of time the surface temperature must he T aw , that is, 

Kt— >co = 1 (5) 

A boundary condition on the temperature gradient at the leading edge 
may he established as follows. Assume that the plate is of unit width. 
Then the frontal area is 1 X h, a finite constant. The condition of no 
heat transfer through the frontal area then requires that the temperature 
gradient there he zero, that is, 



= 0 


Rx-o 


( 6 ) 
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The process of finding a solution of equation (3) proceeds as fol- 
lows. At small values of time a good approximation (zero order) to 
equation (3) is obtained by putting K equal to zero in the second term 
of the right member; thus, 


ac _ a 2 + b 

^ " Sr x 2 b * 1 / 2 


It can be verified that 


(7) 


K 


~ t 5 ^ F(t}) 
fa 


where 


( 8 ) 


T] 


Jk_ 

at 1 / 2 


satisfies equation ( 7 ) provided F 


satisfies 


d% [ i\ dF 3 F _ 1 
dT] 2 2 dt] 4 


(9) 


In order to satisfy the boundary and initial conditions, F must remain 

finite as T) > 00 and F'(o) must equal zero. Either a power-series 

method or numerical integration of equation (9) will give such a solu- 
tion. An iteration process may be used for obtaining higher-order 
approximations to the solution of equation (3); for example, the first- 
order approximation is obtained by substituting the zero-order approxi- 
mation (8) for K in the second term of the right member of equation (3). 
An equivalent process, suggested by the form of the zero-order approxi- 
mation, is to write the solution of equation ( 3 ) in the form 


K = 



t 3n/4 


FnOl) 


( 10 ) 


where 
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substitute equation (10) into equation ( 3 ), collect terms, and equate 
the coefficients of the powers of b/\fa to zero. This solution is pre- 
sented in appendix A, where the result is shown to he a set of ordinary 
differential equations for the functions F. The convergence of equa- 
tion (10) for large values of q is easily demonstrable. The actual 
results for the first three terms indicate that equation (10 ) is con- 
vergent for small values of q. Values of F which permit equation (10) 
to satisfy conditions (^) and (6) may be calculated in a straightforward 
way (see appendix A). The functions F n have the limiting values 


F n (q — *“) 


(-if* 1 

nj l,”/ 2 


and 


F n (q — >0) = C n 


where the coefficients C n are constants positive for odd n and nega- 
tive for even n. 

In the case of zero conduction (a = 0) the solution of equation ( 3 ) 
which satisfies the boundary and initial conditions is 


K_ 



( 11 ) 


Equation (ll) is also the asymptotic solution of equation ( 3 ) with con- 
duction (a ^ 0) in the limit of large q and, therefore, correctly gives 
the temperature at finite values of R x and sufficiently small values of 
time. This result is plausible physically, since at sufficiently small 
values of t ime and finite rate of heat transfer (r x finite), conduction 

effects should have affected the temperatures by a negligible amount. 

The limit of large q also corresponds to large R x and finite time. 
Therefore, equation (ll) indicates that the effect of conduction is 
negligible at sufficiently large values of R x . The functions F-j_, F 2 , 

and F^ as derived in appendix A are plotted in figure 1 against q. 

For comparison, the asymptotic functions (zero conduction) are also 
shown. These results indicate that for values of q larger than about 5 
the F n functions are negligibly different from their asymptotic values 

and hence equation (10) is negligibly different from equation (ll). 
Therefore, for R^. > 5a\/t, the effect of conduction on the temperatures 

is negligible. 
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The solution for small values of tj is perhaps of more interest. 

As T| goes to zero, the values of F n become the constants C n . There- 
fore, the solution (eq. 10) at the leading edge (R x = 0 and q = 0) 
becomes 


where 



( 12 ) 


At small values of time, 
as the 3 A power of the time. 


the temperature of the leading edge varies 
A curve of the form 



0.2041 ~ 
va 



- 9.447 


( 13 ) 


fits the leading-edge solution (12) in such a way that expanding equa- 
tion ( 13 ) allows the first two terms of equation (12) to he given cor- 
rectly and = 1.605 (0.79 percent difference). 

The functions F 2 , and F^ were calculated and are plotted in 

figure 1. In order to present instantaneous temperature distributions 
over the flat plate, it is convenient to define a dimensionless time 
parameter 


a 




and a dimensionless distance parameter 


P 


A/5 


These parameters are related to Tj in the following manner: 

n = 

Figure 2 shows the temperature distribution on the plate at a few chosen 
values of a obtained by using the first three terms of the summation in 
equation (10). The vertical marks on these curves indicate the contribu- 
tion from the third term in the summation and, to some extent, the accu- 
racy. The fourth term in the summation would give a negative contribution. 


P 

a 2 /5 
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The Conical Shell 


The heat conduction equation for the conical shell is derived in 
the following manner. The principle of conservation of heat is applied 
to a ring section of the conical shell of thickness ds parallel to the 
conical surface. Temperature gradients normal to the conical surface are 
assumed to he zero. When the shell is of constant thickness h and the 
axially symmetrical rate of heat input per unit area through the outer 
conical surface is the value of q c given in equation (2), the temper- 
ature of the shell is 


8t 

Vo* * - 


% \Ss 2 s &) 8 v B.V2 


or 


+ j/2L(i-k) 

St VaR a 2 R s fegj Bb 1/2 


(It) 


where s is the distance along the conical surface from the vertex and 


2 = k h 

’ v 2 


h 


0.33 


k g 


U Pr 1 /^ 
v h 


K = 



The fact that the initial temperature is Tj_ and the final temper- 
ature is T aw gives the boundary conditions 

(15) 

(16) 


* 1=0 = 0 

and 

= 1 
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Since the frontal area of the conical shell is zero, the condition 
of zero heat transfer through the frontal surface does not lead directly 
to a boundary condition on the temperature gradient at the nose. How- 
ever, the condition that 

K remains finite for all values of t and R s (17 ) 


(which, incidenta l l y, is implied in equations (15 ) and (1 6)) is adequate 
to select a unique solution. It is interesting to note that the solution 
selected by condition (17) has a zero temperature gradient at the nose. 
The assumption that the conical-shell thickness h is constant all the 
way to the vertex is an obvious mathematical fiction for values of s 
less them h/tan 0 where 0 is the cone half angle. This error in the 
analysis of the conical shell is in such a direction as to give a rate 
of temperature rise that is less than the actual value, that is, too 
large a correction due to the effect of heat conduction. (This fact, 
partly at least, prompted the analysis of the next section. ) 

The analysis of the conical shell (appendix B) is very similar to 
that for the flat plate. The solution is written as 



where 


Rc 


at' 


1/2 


The zero-conduction (or large q) solution is 


-a- 


bt 


Kc = 1 - e 

The nose solution (q = 0) is 


R, 


,V 2 


(19) 




( 20 ) 
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where the coefficients in the first three terms of the summation are 

d-l = 1.15622 
d 2 = -0.69618 

Dj = 0.29080 


The expression 



|l + 0. 04801 



-24.08 


( 21 ) 


fits the nose solution (eq. (20 ) ) almost perfectly for the first three 
terms. The functions G-^, Gg, and were calculated and plotted in 

figure 3. 


In order to present Instantaneous temperature distributions over 
the conical shell, it is convenient to define a dimensionless time 
parameter 


_ \/3 ht 5 ^ 


and a dimensionless distance parameter 


P = 


R^V 5 


related to q hy 


4 = 


P 

a 2 / 5 


Figure 4 shows the temperature distributions obtained over the conical 
shell for a few chosen values of a by using the first three terms of 
the summation in equation (l8). The vertical marks on these curves 
indicate the contributions from the third term. 
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The Slender Solid Cone 

The transient one -dimensional heat-conduction problem applies also 
to the solid cone in the limit of small cone angles. The heat-conduction 
equation is derived by applying the principle of conservation of heat to 
an element of the solid cone, which is a spherical shell with center at 
the cone vertex of radius s and thickness ds. For the axially symmet- 
rical case, q c from equation (2) is the rate of heat input per unit area 

through the conical surface. Temperature gradients tangential to the 
spherical shell element are assumed to be zero. The resulting equation 
for the temperature of the slender solid cone is 


or 


where 


8t f c£r 2 8r\ sin e ^c 
bt b V ds 2 + 8 8s / + 1 - cos 9 s 


8k _ 2 / b^K 2_ cK 

at = a Us 2 + Rs *BJ R e 


t>C 


£- — ]+ (1 - K) 

” ' ' 3/2 


,2 _ 


k b U 2 


P,. c 


b b 


t> c = 0.33\/3 ^ Pr^g 

v 2 


( 22 ) 


sin 0 

g = 

1 - cos 0 


K = 



and 


0 is the half vertex angle of the cone. 
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Since the initial temperature is Tj_ and the final temperature is 
T aw , the solutions of equation (22) must satisfy the conditions 


and 


*t=0 = 0 


K. __v 

t 7 OO 


» 1 


( 23 ) 

( 24 ) 


and the condition that: 


K re main s finite for all values of t and R s 


( 25 ) 


The condition (25 ) specifies a unique solution which has a gradient at 
the nose varying inversely as the square root of R s . 

The solution is straightforward (appendix C), and may he written as 


where 


CO 


n=l 



HnOl) 


n 


R s 



The zero conduction (large T]) solution is 


Koo = 1 



(26) 


(27) 


The nose solution is 


Kno 


se 



En 


(28) 
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where 


E 2 = 1.95940 
Eg = -2.20696 

E^ = 1.91106 

The expression 



1 - 1 + 0. 2935 


, 3/2 


,1/V 


-6.681 


(29) 


fits the first two terms of the summation of equation (28) and gives the 
third term within 2.0 percent. The functions Hj, Hg, and were 

calculated and plotted in figure 5- 

In order to present instantaneous temperature distributions over the 
slender solid cone, it is convenient to define a dimensionless time param- 
eter by 


bet 1 / 4 

a = 

and a dimensionless distance parameter 


(3 = 


Esb, 


B U C 


related to T] by 


1 



Figure 6 shows the temperature distribution over the slender solid cone 
for a few chosen values of a obtained by using the first three terms 
of the summation in equation (26). The vertical marks on the curves 
indicate the contributions from the third term of this summation. 
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DISCISSION OF RESULTS 


Figures 2, k, and 6 present the temperature distributions over the 
flat plate, conical shell, and slender solid cone, respectively, as a 
function of the parameters a and f3. On the basis of expressions (l) 
and (2) for the heat transfer through the laminar boundary layer, a 
and p are known functions of the flow properties of the air and the 
variables of the problem in each of the three cases. However, since in 
each of the cases analyzed the flow is steady just outside the boundary 
layer, the only real assumption is that the temperature' end Reynolds num- 
ber dependence of the heat- transfer expressions is correct. 

Generally it is evident that the conduction effects are important 
near the nose or leading edge. In the most forward region the temper- 
ature is always lower than the zero-conduction temperature. Immediately 
behind this region is a second region, for which the rearward boundary 
is not well defined and where the temperature is higher than the zero- 
conduction temperature. The boundary between these two regions moves 

rearward as time increases. Because of the R^ and R^ factors in the 
mass elements, the temperatures for the conical shell and for the slender 
solid cone are negligibly different in the second region from the zero- 
conduction temperature. 

Nonweiler (ref. 3) calculated the equilibrium temperature of the 
flat plate by taking into account internal conduction on the assumption 
that the heat loss by radiation is of such importance that the equilibrium 
temperature is greatly different from the adiabatic wall temperature. He 
obtained an equilibrium tenperature distribution which is quite similar to 
the instantaneous distributions of figure 2. However, Nonweiler' s result 
derives from a different cause. A conduction effect is present in his 
time -independent case only because the adiabatic wall temperature is never 
reached because of radiation heat transfer, an effect omitted in the 
present analysis. Nevertheless, qualitatively, the effect of heat con- 
duction near the leading edge is s im ilar in the two cases. 

For the flat plate and the conical shell the tenperature at the 
leading edge or nose varies initially as the 3/^ power of the time, with 
a coefficient which is nearly the same for the two. Specifically, 


First Term of Flat-Plate Solution ^1 _ 1 

First Term of Conical-Shell Solution ^3 Dp 1. 04 

If expressions (13 ) and (2l) are assumed to be correct for the nose 
teiperatures, the ratio of the time required for the flat-plate leading 
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edge to attain 1/2 of its total temperature rise (nose half -rise time) 
to the time required for the conical-shell nose to attain l /2 of its 
total rise is 1.08. It must he noted, however, that the assumption of 
a constant thickness. all the way to the vertex of the shell results in 
a calculated rate of nose temperature rise less than the actual value. 

The nose temperature of the slender solid cone varies initi ally as 
the 1 j\ power of the time. Since the solid-cone solution depends on the 
flow and the body properties in a way different from that in the case of 
the flat plate and conical shell, it is not possible to compare them 
except for specific cases. The obvious and striking difference is that 
the temperature gradient is infinite at the nose and that the tempera- 
ture falls off much more rapidly for the solid cone than for the plate 
and shell. 

It is interesting to estimate the extent of the forward region over 
which heat conduction effects are significant. For the flat plate the 
second region (temperature higher than zero -conduction temperature) is 
significant. Approximately, q = 3-5 defines the extent of the con- 
duction region; thus : 

For the flat plate: 

B *x 

”° ' JR = ^ = 3 ' 5 

For the conical shell and the slender solid cone, the second region is 
not significant. The position at which the actual temperature equals 
the zero-conduction temperature is taken as defining the significant 
region. Approximately, these positions hre: 

For the conical shell: 


% - — ^ 7 - * 1-75 

a 2 / 5 at 1 /* 


For the slender solid cone: 


— = is — = 2.0 

° a* at 1 / 2 
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Solving these defining criteria for Reynolds number gives: 
For the flat plate: 


^/2 


= 5 * 5 (^) V Mtl/2 = 9 * 5 X 105 


For the conical shell: 


vV2 




= 1 . 75 /!^-^ - Mt 1 / 2 = 4.7 XlO 5 Mt 1 / 2 

ww v 


For the solid slender cone: 


R = 2. 0 /— ^ ^ - Mt 1 / 2 = 5-3 X 10? Mt 1 / 2 
S \%%) V 


where the final expressions (here and subsequently) are calculated for 
standard air and copper. 

If the approximate expressions (13 ), (21), and ( 29 ) for the forward- 
edge temperatures are used, values of a corresponding to the nose half- 
rise times are: 

Flat plate: 


a = 0.37 = t 5 /^ = 0.000357 — 

Va h 

Conical shell: 

a = 0.6l = ^ t 5 ^ = 0.000613 ^ 

h 


Slender solid cone: 


a = O .37 = 


, 5/2 


t 1 ^ = 0. 0200g'/M 
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Combining these results gives for the Reynolds number extent of 
the significant region at nose half -rise time: 

Flat plate: 



aV5 ,(0.57) 2/5 


,2/5 


= Pr 


-2/9/057 c 

\0.53 kg v 



9.66 X 10 ^ M 2 /^ 2 / 5 


Conical shell: 


R 


>/3 


2/3 


-(0 


. 61 ) 2/3 - Pr -2/9/061 C J 69 x 10 7 ^2/5 


S ° >#3 \ 0 - 33 # \ ^ 7 


Slender solid cone: 


R s. - — gf 0 - 37 ) 2 
0 b c 2 


0.37 k b\ 1 2A2 X 10 


8 


^033^3 k g/ ^2/3^ g ‘ 


The nose half -rise times are: 

Flat plate: 

M a 2 /^ 


t Q = (0.37)' ^ — 7 - = 1-99 x 10^ - 


k /3 


>/3 


M 2 /3 


Conical shell: 


(0.6l) 4 ^ a 2 /? 

5 2/3 b ty3 


1.01 X 10 


,4 & 


4/3 


m 2 /3 


Slender solid cone: 


t 0 - (0.37) 


a c 


2.07 x 


10 5 


s 


V 
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It is evident that the results depend upon the properties of the 
material of which the bodies are made primarily through the conduction 
parameter a. The heat-transfer parameters h and h c depend on the 

product p^c^, varies only slightly from one metal to another. 


However, the parameter a is proportional to 



For the flat 


plate and the conical shell the nose half -rise time varies as the 
l/3 power of and the extent of the conduction region varies as the 

2/3 power of k^. On the other hand, for the slender solid cone the 
nose half -rise time varies as the cube of k^ and the extent of the 


conduction region at nose half -rise time varies as the square of k^. 


A serious limitation of the present analysis results from the 
assumption of uniform temperatures in the bodies in a direction perpen- 
dicular to the stream direction, especially in the case of the slender 
solid cone. Undoubtedly the solid-cone results are valid only in cases 
of rather small half- vertex angles. The geometry factor g varies, for 
sufficiently small vertex angles, as 2/8. Because of the strong depend- 
ence on the factor g, with cone angles for which the analysis is reason- 
ably correct, equilibrium conditions near the vertex are attained in 
extremely short time intervals. The error is not so serious for the 
flat plate and the conical shell, although as the thickness h increases, 
the error increases. 

For the flat plate and the conical shell the extent of the signifi- 
cant region and the nose half -rise time vary only slightly with Mach 
number and thickness. The dependence of the results upon Mach number 
and kinematic viscosity always occurs in the combination cM/v, so that 
at constant temperature similar situations result if the Mach number is 
proportional to the kinematic viscosity. 

Ho suitable experimental data are available to check the results of 
this analysis. 


COHCHJSIONS 


The analysis of the transient heat- conduction effects for the semi- 
infinite bodies the insulated flat plate, the conical shell, and the 
slender solid cone, subjected to steady flight conditions beginning at 
zero time and with the only heat transfer that through a laminar boundary 
layer, shows that the effect is important only over the most forward por- 
tions of the body. For all three bodies, the extent of the conduction 
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region Increases with the square root of time and linearly with the Ife.ch 
number. All results depend upon Mach number and free -stream kinematic 
viscosity as the quotient of the two, so that similar results are obtained 
when the Mach number varies directly as the kinematic viscosity. 

For the flat plate and conical shell in typical flight conditions, 
the time required for the nose to attain one half the total temperature 
rise (nose half -rise time) is of the order of a fraction of a minute and 
varies inversely as the 2/3 power of the Mach number. The extent of the 
conduction region at nose half -rise time is of the order of inches and 
varies linearly with the Mach number. 

For the slender solid cone the nose half -rise time and extent of 
conduction region at nose half -rise time depend strongly on the cone 
half angle. With cones sufficiently slender that the analysis is reason- 
ably correct and in the same flight conditions, the nose half -rise time 
and extent of the conduction region are one (or more) order of magnitude 
less than for the plate and shell. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 13, 1955* 



2k 


NACA TN 3058 


APPENDIX A 


SOLUTION OF FLAT- PLATE PROBLEM 


The solution of equation (3) 



(Al) 


(A2) 


where 


T) = 



Substituting equation (A2) into (Al) and equating coefficients of 
the powers of b/i/a yields 


V + i 


V + i F 2 ' 


3 p _ 1 
ir F i-~= 




3 F = £1 
2 




(A3) 


J n-1 


1 11 + 1 p ' _ 3n 
n 2 ^ *n ^ 



4T 
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The asymptotic (large T|) solutions for the F functions are 



(-d m1 


which when substituted into equation (A2) give the asymptotic solution 
(zero-conduction solution): 


Koo = 



m 


The ordinary differential equations for the F functions may be solved 
in sequence; F^ may be written as the sum of a particular solution and 

a power-series homogeneous solution 


Ft = - — T] 5 / 2 + C-i 


1 + \ r\ 2 + . . . + 

8 128 


(-l) r+1 d( ljljr)^ + . . 


2r 


(A5) 


where the recursion formula for the coefficients is 

d(l,l}i*l) = ■ |4r ~ d(l,l;r) 

4(2r + l)(2r + 2) 

The quantity F^'(o) = 0 which is required by the boundary condition (6). 
Proceeding with similar solutions for F2 and F^ yields 


Fo = - 


- n 5 + C x 
9 1 


i + 


d(2,l;l)Tj7/2 _ . 


.(-l) r+1 d(2 ; l;r)Ti"r + 


1 + d(2,2;l)T] 2 + d(2^2;2)t)^ - 


d(2,2;5)tl 6 + . - .(-l) r d(2,2;r)ri 2r + . 


(A 6) 



where 


d(2, ljr+l) 


4d(l,ljr+l) + |4r - 3 |d(2,ljr ) 
(4r + 7 ) (tar + 5) 


d(2,2;r+l) 


\_2r - 3] 

2(2r + l)(2r + 2) 


d(2,2;r) 


and 


F-r = - 


8 

567 


^ + c n 


l 3 + d(3,lj l)r) 3 - 


d(3,l;2)ri7 4- . . .(-l) I ^ :i nd(3,ljr)Ti 2r+ 5 + . . . 


| T] 3 / 2 + d(3,2jl)T] T / 2 + d(3,2}2)ri 11 / 2 - 
3 

/ _ 4r+3 

a(3,2}3)Ti 1 5/ 2 + . . .(-l) r d(3,2jr)n 2 + . . . 


1 + d(3,3;l)n 2 + a( 3 , 3 ; 2 )T] 4 + d(3,3>3 )vfi 


a(3,3j^h 8 + . . .(-l) r+1 d(3,3;r)Ti 2r + . . . 


where 


d(3,l;iH-l) = 

d(3,2;r+l) = 
d(3,3>r+l) = 


4d(2,l;r+l) + \kr - 3jd(3,l;r) 
k(2r + 5)(2r + 4) 

l)-d(2,2;r+l) + |4r - 6|d(3,2;r) 
(4r + 7)(^r + 5) 

| 4 r - 9 f 


k(2r + l)(2r + 2) 


a(3,3;r) 
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The series expansions for the F's (that is, equations (A 5 ), (A6), 
and (A7)) may be used to evaluate the F's, or the differential equa- 
tions (A3) may he integrated numerically. The series expansions were 
used for the calculation presented herein. For n > 3 "the series 
method probably involves more labor but can yield greater accuracy. In 
either case, the boundary condition (6) requires that 

V(o) = 0 


and the constants C n 


F n (0) = C n 


must be determined so that each F has its proper asymptotic (large t|) 
value. 


*n(l 


n.'T] 11 /^ 


Each new F involves one new C. Utilizing the fact that at suf- 
ficiently small values of time the solution (eq. (A2)) reduces to one 
term permits Cp to be calculated analytically. 


The method consists in constructing a solution of the one- 
dimensional heat-conduction problem for boundary conditions consistent 
with the present problem at small values of time and then evaluating this 
solution at the leading edge of the plate. A well-known solution of 


Is 


dT 2 8 2 T 

— = a 

* 8x 2 


(A8) 


T*(x,t) 


1 e -(x-£) 2 Aa 2 t 

2a^/rtt 


(A9) 


This solution corresponds to an initial temperature distribution which 
is zero everywhere except at x = | and is infinite at x = | in such 
a way that 

n+» 


-00 


T*(x, 0)dx = 1 
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In other words, equation (A 9 ) corresponds to an initial temperature 
distribution given by the Dirac 5 function, b(x - g). The essenoe 
of the method is to interpret equation (A9) as the solution corresponding 
to injecting a unit quantity of heat at x = | at t = 0 . 

The construction of a solution satisfying the conditions of the 
present problem (at small time) proceeds as follows. In order to insure 

8t 

w- =0 at x = 0 , the input is made symmetrical about x = 0 : 


T(x, g,t) 


1 

2a\/jrt 


g -(x-|) 2 Aa 2 t + e -(x+|) 2 Aa 2 t 


(A10) 


Equation (AlO) corresponds to unit heat input at t = 0 at x = g 
and x = -g. Then 


T(x,t) 



poo 

/ L 

e -(x-g) 2 Aa 2 t + e -(x+g) 2 Aa 2 t 

Jo 'lift 



d| 


(All) 


corresponds to putting 



units of heat per unit time at 


t = 0 


into the interval dg at g. Then, if the input rate is assumed constant 
in time for sufficiently small times, and is integrated over a small time 


t 


F(x,t) = 


It poo 


2aft Jo J 0 \/il/tL 


■(x-g) 2 Aa 2 t ^ -(x+g) 2 Aa 2 t 


+ e 


dg dt 


(A12 ) 


the expression for the temperature at x and at a small time t is 
obtained. At the leading edge (x = 0), equation (A12) reduces to an 
expression which can be evaluated 


F(0,t) 


* rr 2 e -g 2 Aa2t 
2avGr Jo J 0 \Uft 


dg dt 


Integrating first with respect to g gives 



2 r(A 


(A13) 


(A14) 


F(0,t) = 


dt 
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and., finally, integrating with respect to t yields 


F(0,t) = t 3 ^ 

Va 


^ r (t) 

3 ^ 


(A15) 


Comparing equation (AI 5 ) with the first term of equation (A2) at x = 0 
(tj = 0) yields 



1- 9285 


The coefficient Cq_ could also he determined by equating the heat 
inputs in the cases of conduction and zero conduction (at small times ) 
over a sufficiently large forward portion of the plate. This procedure 
would result in equating the integral of equation (A5 ) at, say, tj = 6. 0 
to the integral of the asymptotic (zero -conduction) value of F^ the 

same tj. The value of C is checked by F-j_ approaching its proper 
asymptotic value. The C n 's for n > 1 are determined by choosing 
them so that the corresponding F n *s have their proper asymptotic 
values. 
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APPENDIX B 

SOLUTION OF THE CONICAL-SHELL PROBLEM 
Substituting equation (l8) 



atV2 


into equation (lb) 

+ Ja.(i.K) 

\ 8r s 2 R s / Rg 1 / 2 

and equating the coefficients of the powers of b/\/a yields 



(Bl) 


(B2) 


(B3) 


G n " 



3n = ^-1 

" 4 ' fn 


The asymptotic (large q ) solutions for the G functions are 




0 = 


(-l ) n+1 

n.' r) n / 2 


(B4 ) 
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(the same, incidentally, as for the F functions). 'When the asymptotic 
values of the G functions (eq. (b 4)) are substituted into the solu- 
tion (Bl), the asymptotic (zero -conduct ion) solution (eg. (19)) results. 

The Boundary condition (17) (K remains finite for all values of t 
and R s ) requires that 

G n remains finite for all values of T| and n (55) 


Series solutions for the G functions satisfying condition (B5) are 


°i 




+ Di 


1 + 3 2 3 4 + 

16 102k 


(-l) r+1 d(l,l;r)ri 2r + 


(B 6) 


where 


d(l,l}r+l) = 3 L d(l,l;r ) 

i|-(2r + 2 r 


Similarly, 


G 2 = ^ + D ! | 1 1 3/2 + d(2,l> l)i)7/2 _ 


4r+3 


d(2,l}2)ii 11 / 2 + . . . (-l) r+1 d(2,ljr)T] 2 


Dr 


1 + d(2,2jl)T] 2 + d(2,2;2)^ - 


d(2,2;3)Tl 6 + . . .(-l) r d(2,2jr)T 1 2r + 


where 


d(2,l;r) = ^ d(l ' 1?r) + - 7|a(2,l;r-l) 

(4r + 3) 2 

d(2,2;r) = l 21 ’ - " 5 . 1 . d(2,2;r-l) 

Qr 2 
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^9/2 + d — il 5 + d(3,ljlh 5 - 
6561 |_ 8 l 


4(3,1j2)t|7 + 


(-l) r+1 d(3,ljr)Tl 2r+5 + . . . + 


T>2 [I I! 5 / 2 + d(3,2jl)ri^/ 2 + d(3,2}2h rL / 2 - 

4r+3 —1 

d(3,2}3)il 15 / 2 + • • .(-l) r d(3,2;r)Tj 2 + . . . 

D 3 Jl + d(3,3;lh 2 + d(3,3j2h 4 + d(3,3}3h 6 - 
d(3,3}^)Tl 8 + • - .(-l) r+ 1 d(3,3jr)Ti 2r + . . . 


d(3,i;w-i) 


d( 3 , 2 jrt-l) = 


4d(2,i;«-i) + l 4 ** - 3la(3A;r) 

k(2r +5 ) 2 

lfd( 2 , 2 $iH-l) + |4r - 6 |d( 3 , 2 ;r) 

(4r + 7) 2 


d(3,3ji+l) - |4r - 9l ^ d(3,3}r) 

4(2r + 2 ) 


No straightforward analytic method for determining the constant 

was devised. The D coefficients were determined "by choosing them so 
that the series expressions for the values of Gq were close to their 

asymptotic value at t\ = k. 5 . These constants appear in the nose solu- 
tion (x = 0 ) since (3^(0) = E^. 
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APPENDIX C 


SOLUTION OF THE SLENDER-SOLID-CONE PROBLEM 


Substituting equation (2 6 ) 


-S 


where 


n = 


Rp 


at x /2 


into equation (22) 




:(l - K) 


\ ^ 2 ■ R s Sr r ) v 3/2 

and equating the coefficients of the powers of b c /a3/2 yields 




v + (f*0%’ -1*2 = J75 




V + (f + iH' ' f * 


Hn-1 

^5/2 


(Cl) 


(C2) 


(C3) 


The asymptotic solutions for the H functions are 


n+1 


HnC 1 ! - = 


(-D 

n.'T|3 n / 2 


(C4) 
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When the asymptotic values of the H's (eq. (C4)) are substituted 
into the solution (Cl), the asymptotic (zero-conduction) solution (27) 
results . 

The boundary condition (25 ) (K remains finite for all values of t 
and R s ) requires that none of the H functions diverge as Tj — > 0. 

The series solution for H]|_ satisfying this boundary condition is 


H 1 = + E 1 


d(i,i } 2)n 4 + 


1 4- d(l,l)l)7] 2 - 

. . ,(-l) r+1 d(l,l;r)T] 2r + . . . 


where 


d(l,l;r+l) = 


|4r - l| 

4(2r + 2)(2r + 3) 


d(l,l;r) 


The series solution for is 


H2 = - — i\ + 


| R 1 / 2 + d(2,ljl)T]5/2 


+ . . .(-l) r+1 d(2,ljr 2 + . 


+ 


1 + d(2,2;l)r] 2 - d(2 ; 2>2)T] 2f + . . . (-l) r+1 d(2,2;r )ii 2r + 


where 


d(2,l;r + l) = ^dA?r+l) + [hr - lfd(2,l?r) 

(hr + 5)(4r + 7) 


d(2,2jr+l) 


|2r - 1[ 

2(2r + 2)(2r + 3) 


d(2,2;r) 
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The result for is 




E-, 


| 1 \ + d(3Ajlh 5 

3 


d(3,lj2)T]5 + . . .(-l) r+1 d(3,l}r)T] 2r+1 + . . . 


2 r+l 


E r 


| T) 1 / 2 + d(3,2jl)T]5/2 


4r+l 


d(3,2j2)r}9/2 + . . .(-l) r+1 d(3,2;r)T 1 2 


+ • • 


=5 


1 + d(3,3il)n 2 - d(3,3}2)n 4 + . . .(-l) r+ 1 d(3,3ir)T 1 2r + 


where 


d( 3 ,l;r+l) = 


4d(2,l>r+l) + |4r - l|d( 3 ,l;r) 
4(2r + 3 )( 2 r + 4) 


d( 3 , 2 jr+l) = Ud ( 2 , 2 ;r+ 1 ) + j 4r - g|d( 3 , 2 ;r) 

(4r + 5 )(4r + 7 ) 


d(3,3>r+l) = — - |4r ' 31 d(3,3jr) 

4(2r + 2 )( 2 r + 3) 


As T] — > 0, the H functions take the values ^ which, 

when substituted into the solution ( 26 ), gives the nose solution ( 28 ). 
The constants E^ (up to Ej) were determined by choosing them so that 

the series expressions for the H functions were close to their asymp- 
totic values at rj = 4 . 5 . 
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Figure 2.- Temperature distributions over the flat plate at three values 



of the time parameter a. a 
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Figure 5.- The first three functions in the slender- solid- cone solution 
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